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We engineer the fast rotation of an effectively one-dimensional ion trap for a predetermined 
rotation angle and time, avoiding the hnal excitation of the trapped ion. Different schemes are 
proposed with different speed limits that depend on the control capabilities. We also make use of 
trap rotations to create squeezed states without manipulating the trap frequencies. 

PACS numbers: 


I. INTRODUCTION 

A major challenge in modern atomic physics is to develop a scalable architecture for quantum information processing. 
A proposed scheme to achieve scalability is based on shuttling individual ions or small groups of ions among processing 
or storing sites Apart from shuttlin g [3 -[lTI|. other manipulations of the ion motion would be needed, such as 

expansions or compressions of the trap [l^ Il3l|. separating or merging ion chains [iM| , and rotations [I3- All these 
basic dynamical operations should fulfill two seemingly contradictory requirements: they should be fast, but free from 
final motional excitation. Shortcuts to adiabaticity for “fast and safe driving” have be en p roposed for several of these 
elementary operations M,[il and have also been implemented experimentally [islIl9j| . 

In this paper we study rotations of a single ion as depicted in Fig. [1] Our aim is to inverse engineer the time- 
dependence of the control parameter(s) to implement a fast process, free from hnal excitations. We assume for 
simplicity that the ion is trapped in a linear, harmonic trap, tightly conhned in a radial direction so that it moves 
effectively along a one-dimensional axial direction, hereafter “the line”. The trapping line is set horizontally and is 
rotated in a time tf up to an established hnal angle = 7r/2 in all examples) with respect to a vertical axis that 
crosses the center of the trap. Such an operation would be useful to drive atoms through corners and junctions in a 
scalab le q uantum processor [laiiii. It is also a hrst step towards the more complicated problem of rotating an ion 
chain |l7l [^. [^. which would facilitate scalability in linear segmented traps, and be useful to rearrange the ions, 
e.g. to locate a cooling ion at the right position in the chain [i3- 


t.f 



FIG. 1: (Color online) Schematic representation of the rotation process. The ion is confined along a line (where it is subjected 
to an effective one-dimensional -longitudinal- potential), which is rotated by an angle 0 up to 0/ in a time tf, so that the final 
state is not excited. 


Electronic address: mikel.palmero@ehu.eus 







2 


We shall first find the classical Hamiltonian. Let s denote a point on the line, s may take positive and negative 
values. A time dependent trajectory s(t) has Cartesian, laboratory frame components x = x(s,t), y = y{s,t), 

X = scos(0), y = ssin(0), (1) 

where 0 = 9{t) is the rotation angle. The kinetic energy is K = ^m(x^ + where m is the ion mass, and the 
potential energy is assumed by now to be harmonic, (this will be relaxed below and in Sec. |n|, where wg is the 

angular frequency of the external confining trap in the (longitudinal) direction of the line. This gives the Lagrangian 

L = 2^^^ ~ 

up- = Pl-P. (3) 

Note that the angular velocity of the rotation 0 must be real but could be negative, whereas lP may be positive or 
negative, making uj purely imaginary in the later case. Unless stated otherwise, the following physically motivated 
boundary conditions are also assumed: the initial and final trap should be at rest, and we also impose continuity of 
the angular velocity, 


0(0) = 0, 0(t/) = 0f, (4) 

0(0) = 0{tf) = 0, (5) 

a;(0) = Lo{tf)=uJo, (6) 

where the last line follows from the second one using Eq. By a Legendre transformation we finally get the 

Hamiltonian^ 

TT • 1 ,2 1 2 2 

H = s— - L = -ms H —mPs . (7) 

ds 2 2 ^ ^ 

At this point, we quantize this Hamiltonian by substituting ms by the momentum operator p and by considering s 
as the position operator, which becomes a c-number in coordinate representation. 


u ^ 2 P 2 2 

H = -—p + —moj s . 
2m 2 


( 8 ) 


We will from now on work with this quantum Hamiltonian (possibly with a more general potential) and corresponding 
quantum states. It represents formally a harmonic oscillator with time-dependent frequency, but there are significant 
differences with an actual harmonic oscillator when the inverse engineering of uj(t) is considered. For an actual 
harmonic oscillator a fast and safe expansion or compression in a time t / should take the system from an initial value 
to a final value of ui without final excitation, in principle without further conditions. By contrast, in the rotation 
process, according to Eq. the initial and final effective frequencies are the same, but the conditions in Eqs. 
and ([5]) must be satisfied. This implies an integral constraint on w. 


0{tf) = / 0dt' 

Jo 



(9) 


where the square root branch should be chosen to satisfy continuity. One further difference is that in a physical 
expansion/compression uj{t) is controlled directly whereas in the rotation there are several options. If Wg is constant, 
only 0{t) is controlled, so that w(t) is an ‘effective’ frequency. In general both wg and 0 could be controlled as time- 
dependent functions, see the next section. As for the final excitation, the expression for the energy of a state that 
begins in the n-th eigenstate of the trap at rest can be found making use of the Lewis-Riesenfeld invariants [l3, , 

see the corresponding time-dependent wave function in the Appendix, 




{2n + l)h 
4a;(0) 


- P{t)P 



( 10 ) 


^ This is easily generalized for a potential U{s), not necessarily harmonic, as H = U{s) — 
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Here 6 is a scaling factor, proportional to the width of the invariant eigenstates, that satisfies the Ermakov equation 

b + ( 11 ) 

To avoid any final excitation, it is required that 

b{tf) = l, b{tf)=0 (12) 

for the initial conditions 6(0) = 1, 6(0) = 0. The boundary conditions for 6 and Eqs. (I4I5I6I) imply that H{0) = H{tf) 
commutes with the corresponding Lewis-Riesenfeld invariant [2^, so that the n-th initial eigenstate is dynamically 
mapped onto itself (but rotated) at time tf. In Eqs. (fTUl) and (fTTl) both the excitation energy and the wave packet 
width are mass independent, so that inverse-engineered rotation protocols will be independent of the species. In the 
following sections we shall analyze different methods to perform the rotation without hnal excitation. 


II. CONTROL OF TRAP FREQUENCY AND ANGULAR VELOCITY 

If both the trap angular frequency ojq and the angular velocity 0 are controllable functions of time, a simple family 
of solutions to the inverse problem is found by setting a 0(t) that satisfies Eqs. (|3]) and ([5]), and compensating the 
time dependence of 0^ with a corresponding change in Wq)^), so that uj^(t) = uj^(O) remains constant during the whole 
process. From the point of view of the effective harmonic oscillator ‘nothing happens’ throughout the rotation, so 
that the effective state remains unexcited at all times. 

We may apply the Lewis-Leach theory of quadratic in momentum invariants to extend the above results to 

arbitrary potentials^. The family of Hamiltonians 

where U is an arbitrary function, and fi depends on time, has the invariant 

where tt = bp — mbs, and Qq is a constant, provided the Ermakov equation 

b + nH=^ (15) 

is satisfied. Consider the simple case ilo = 0, from Eq. m, 

= (16) 

If we set b{t) = 1 as a constant for all times, it follows that H{t) = 0. However, as we saw in the previous section, the 
rotation of a line with the potential U{s) produces in the line frame a centrifugal term —0^s^m/2. To cancel the total 
harmonic term, we have to add to the trap potential a compensating harmonic term, a;gS^m/2, such that = 0^. 
In other words, — 0^ = 0. The resulting Hamiltonian and invariant (in this case they are equal) are simply 

H = I=^ + U{s), (17) 

i.e., time independent. No excitation occurs at any time in spite of the fact that a rotation is taking place. 


^ The theory was first formulated for classical systems in but is applicable to quantum systems as well [^ . Incidentally this means 
that the rotation protocols designed in this paper -in this and the following sections- are valid for classical particles as well. The 
difference appears only when considering which states are valid or not for classical and quantum particles, e.g., when using phase-space 
formulations of quantum states and classical ensembles. 
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For some applications it may be interesting to consider in Eq. (1131) the more general case in which b depends on 
time (for example to achieve a squeezed state), and = uj^ — 0^, corresponding to an auxiliary harmonic term and 
the centrifugal term. The inverse engineering in this case proceeds by designing 9{t), so that 0(0) = 9{tf) = 0, and 
then b(t) obeying the boundary conditions 

6 ( 0 ) = b{tf) = 1, (18) 

6(0) = bitf) = 0, (19) 

6 ( 0 ) = b{tf) = 0 , ( 20 ) 

(or more generally b{tf) = 7 ) that guarantee the commutation between invariant and Hamiltonian at boundary times. 
Once 9 and 6 are set we design the auxiliary harmonic term considering, as before, Oq = 0 in Eq. disi), 

col = n'^+ 0^ =-^+9'^. ( 21 ) 

The auxiliary harmonic term vanishes at both boundary times according to the boundary conditions imposed on 6 and 
0. In fact vanishes as well at the boundary times so that before and after the rotation the atom is confined only 
in the potential U{s). This type of protocols, where both the rotation speed and the potential have to be controlled 
(the latter in space and time) may be quite demanding experimentally. In the rest of the paper we shall assume 
the simpler scenario in which only the rotation speed 9 is controlled, and the trap potential is purely harmonic with 
constant angular frequency luq. 


III. BANG-BANG 

It is possible to perform rotations without final excitation satisfying Eqs. m and dSI) keeping 9 constant or piecewise 
constant. Here we consider the simplest one-step case, 

f 0, t < 0 

9{t) = < c, 0<t<tf . (22) 

[ 0 , t>tf 

Note that Eqs. and (|ni) are only satisfied now as one-sided limits. A bang-bang approach may admitedly be difficult 
to implement because of the sharp changes involved, but it sets a useful, simple reference for orders of magnitude 
estimates of rotation speeds which may be compared to smoother approaches that will be presented later. Integrating 
9 we find 


9f = ctf. (23) 

For a constant 9 = c, u) remains constant from t = 0 to t = tf, and equal to wi = (wq — whereas a; = wq in the 

initial and final time regions. For this configuration, and 0 < t < tf, 


b{t) 

b{t) 


iLOn-OJt 


sin^(a;it) -I- I, 


sin(wi<) cos(wit)(a;o — u;f) 
u}ib{t) 


(24) 

(25) 


to satisfy the boundary conditions 6(0) = 1, 6(0) = 0. The shortest final time to satisfy the conditions (I12|) at t/ is 
tt/wi. From Eq. (1^ this gives the value of c needed, 


whereas 



tf 


f 


TT 

Wl 


TT 



Wo 



(26) 


(27) 


(28) 
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As c < Wo the effect of this bang-bang protocol is to expand the effective trap during the rotation time interval, b 
increases first and then decreases during half an oscillation period of the effective trap. This does not in general 
coincide with half oscillation period of the actual non-rotating trap tt/wq because of the / factor, but it is not too 
different for relevant values of Of. In particular, for Of = 7 r/ 2 , / = 1.118. The maximum of h{t) at tf/2 is precisely 
/. For example, for a frequency wo/(27r) = 2 MHz, this implies a final time t/ = 0.28 fLS. 


IV. OPTIMAL CONTROL BY PONTRYAGIN’S MAXIMUM PRINCIPLE 


While the previous bang-bang method with just one time segment provides a simple guidance, we are also interested 
in knowing the absolute time minimum that could in principle be achieved (even if the “optimal” protocol ends up 
being hardly realizable). Unlike ordinary expansions/compressions, the shortest time protocol for bounded control is 
not of a bang-bang form. To find it we first rescale the time with wq by setting cr = uj^t for t € [0,t/]. Now we set 
the variables 

Xi[(j) 

X2[cf) 

2^3 (O-) 

where u{a) = u{uJot) = with a G [0,wof/]. 

equation (USD and the constraints in (jSD, ([SD and (jlD, and formulate the time-optimal control (OC) problem for 
rotation of a quantum particle on a line as 




Wo VWq 


(29) 


= / u(T)dT, 

Jo 


Then, we can write a control system describing the Ermakov 



such that x\ = X 2 , 

Xn = ^ + {u^ — l)a;i, 
x'l 

x's = u, 

where T = wot/ and the prime is a derivative with respect to cr, with the boundary conditions 

xi(0) = 1, xi(T) = l, 

2 ^ 2 ( 0 ) = 0 , X2{T) = 0 , 

2^3(0) = 0, X3{T) = 0f. 


(30) 


(31) 


Note that we assume that the boundary conditions for it at t = 0 and t = tf can be fulfilled by the use of a sudden 
switch. 


1. Unbounded Control 


We apply the Pontrygin’s maximum principle [2^ to solve the time-optimal control problem (I30L where the 
Hamiltonian is given by 


H[t, X, It, A) — Aq -l- ^ 1 X 2 + A2 — 2 -l- (it^ — 1 ) 2 ;! -I- Aoit, 

Ixi J 

in which A = (Aq, Ai, A 2 , A 3 ) and Aq is either 0 or 1. The necessary condition ^ = 0 gives 

* _ A 3 

2 A22 :i ’ 


(32) 


( 33 ) 
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FIG. 2: 


(Color online) Time evolution 9{t) for the optimal unbounded (a) and bounded (b) control. The rotation angle is 


which minimizes the Hamiltonian and where the co-states Ai, A 2 , A 3 : [0,T] —>■ M satisfy A' 


A'l 

A'2 



—Ai, 

0 . 



dH 
dxi ’ 


i = 1,2,3, i.e., 


(34) 


Solutions are found by solving the equation system composed by Eqs. (l30l) . (l33l) . and (l34l) with the boundary conditions 
at cr = 0 in Eq. m- We have the freedom of choosing different initial values for the Ai(0) to satisfy the boundary 
conditions at T in Eq. (ED). We apply a shooting method and numerically minimize [xi{T) — l]'^+X 2 {T)‘^ + [x 3 {T) — 9f]‘^ 
for these parameters using MATLAB’s ‘fminsearch’ function with 9f = 7 r /2 = 1.5708. The best results obtained are 
for T = 2.2825, which, for the external trap frequency wo/(27r) = 2 MHz used in other examples, implies a final time 
tf = 0.18 fis. The solution found is not exact, {xi{T),X 2 {T),X 3 {T)) = (1.0765,0.0842,1.5650), which might be an 
indication that the system is not controllable. Eigured] (a) shows the time evolution of u for this case following Eq. 
(l33)l but forcing it to be 0 in the boundary times. 


2. Bounded Control 


Now, consider a bounded control with u{a) € [0,1] for all a £ [0,T]. Because the Hamiltonian (15^ is quadratic in 
M, the optimal control that minimizes H is of the form 


= min 


max 


A3 

2A2a;i 



(35) 


The bounded time-optimal control and the resulting optimal trajectory are illustrated in Figure [2] (b). The minimum 
(dimensionless) time that completes the desired rotation is T = 11.9984 and the calculated final state following the 
optimal control is {xi{T),X 2 {T),X 3 {T)) = (1.0083,0.0382,1.5708). For ujo/{2tt) = 2 MHz, the minimal time is 0.95 
/rs. Since u{a) G [0,1], V cr £ [0,r], from (l30l) we see that 0 > 0, and hence the rotation is always forward. In this 
case, X 3 reaches the desired 9f = 7 r /2 at cr = 11.9028, and the control is turned off. Then, the states Xi and X 2 are 
oscillating to reach the desired terminal state (1,0). Figure [5] (b) shows the time evolution of u for this solution. 


V. SMOOTH INVERSE ENGINEERING 

An alternative inversion route that provides smooth solutions is depicted in the following scheme 

9 - —^ 9 -^ a;-^ E[b(tf), b(tf)] 

' • • . rninimize E . ■ ' ‘ 



























7 



FIG. 3: (Color online) Values of the optimizing parameters 04 (thick blue line) and 05 (dashed red line) for different rotation 
times tf. The trap frequency is u}q/{2tv) = 2 MHz, and the final angle 9f — 



FIG. 4: (Color online) Final excitation energy vs final time for the optimized protocol without (solid blue line) and with final 
squeezing ( 7 ^ = 3, dashed red line). The trap frequency is ujq/{2-k) = 2 MHz, and the final rotation angle Of = The initial 
state is the ground state of the trap. 


First, 9{t) is designed to satisfy Eq. (|3|) and Eq. ([5|) with some free parameters. The corresponding 9 and final energy 
are calculated, and the parameters are changed until the minimum energy (and excitation) is found. 

A convenient choice for 0 is a fifth order polynomial ansatz 9 = J2n=o order to satisfy the boundary 

conditions in Eqs. o and we need to fix parameters ao -3 = (0, 0 ,04 + 205 + 39f, — 204 — 805 — 29f). The other two 
parameters, 04 , 05 , are left free in order to satisfy the remaining two boundary conditions in Eq. © and suppress 
the final excitation energy. In practice we solve numerically Eq. m to find the final energy for each pair 04 , 05 , 
and use MATLAB’s ‘fminsearch’ function to find the values of the free parameters that minimize the final excitation 
energy. 

In Fig. [3] the values of the free parameters that result from this process are given, and in Fig. |4]we depict the 
corresponding excess energy with respect to the ideal target state (as in previous examples, a;o/(27r) = 2 MHz). 
Vanishing residual excitations are found for times shorter than half an oscillation period up to a time t/ ~ 0.23 ps, 
not much larger than the unbounded-optimal-control minimum of 0.18 fis. Fig. [5] depicts the difference between the 
ideal value of b{tf) and the actual value, and makes evident the sharp change that marks the shortest time for which 
a solution exists. Since we have limited the possible solutions by imposing a functional form of the function 9(t), 



bang-bang OC(unbounded) OG(bounded) inverse engineering 

tf (ps) 

0.28 0.18 0.95 0.23 


TABLE I: Minimal rotation times for the different methods. Trap frequency tJo/(27r) = 2 MHz. In bounded OC, 0 < 0 < u>o- 











1.2 



0.1 0.15 0.2 0.25 


tf (IJS) 


FIG. 5: (Color online) Difference between ideal and actual value of b at the end of the rotation vs final time for the optimized 
inverse-engineered protocol for rotations without (solid blue line) and with final squeezing ( 7 ^ = 3, dashed red line). The trap 
frequency is ujo/{2tt) = 2 MHz, and the final rotation angle 6f = ^. 



FIG. 6 : (Color online) Values of the optimizing parameters 04 (thick blue line) and as (dashed red line) for different rotation 
times to generate a squeezed vaccum state with 7 ^ = 3. The trap frequency is cuo/(27r) = 2 MHz, and the rotation angle 


this time is larger than the one found via OC. Note also that the shortest final time is slightly better than the one 
provided by the simple bang-bang protocol. Table I summarizes the results. 


VI. WAVE PACKET SQUEEZING 

Consider now a trap rotation with constant trap frequency wq satisfying the conditions BMl) . and b satisfying 

6 ( 0 ) = 1 , 6 ( 0 ) = 0 , 

b{tf) = 7 , 6 ( 1 /) = 0. (36) 

Unlike the previous sections, 6 ends in a value 7 different from 1. 

According to Eq. (IA3I1 . each initial state (/)n(0) will evolve into at t/, where g = g(tf) = 

dt'/b^(t')^ and (j)n,sq is the normalized eigenstate for the trap with angular frequency uj^q = . (This is a 

virtual trap, let us recall that the actual trap has angular frequency wq.) 

A coherent state at time t = 0, 


|a) =e-l“l i^^(o)) 

n V 
n—0 ^ 


(37) 
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will thus evolve into 


OO ~n 

mf)) = e—^ (38) 

n=0 

where a = This is a coherent state for the virtual frequency LOsq and therefore a minimum-uncertainty- 

product state. However, since the actual trap has frequency Wq, it is also a squeezed coherent state with respect to 
the actual trap, |[r, 5]), see [2^, where r = — In 7 , up to a global phase factor. The final and initial coordinate and 
momentum widths are related by = Ap.o/ 7 - We may rewrite the state at time t/ in terms of the 

squeezing and displacement operators as 

= e-“«3/25(r)i7(d)|0) = 5 ])^ ( 39 ) 

where S{r) = \ a and are annihilation and creator operators for the wo-harmonic trap, and D{z) = 

e^a “ is the displacement operator. Note that the phase at t/, arg(Q;), is controllable by means of the g-function 
that depends on the process history, whereas the squeezing parameter I /7 is controlled by the imposed boundary 
condition. If necessary, a controlled tilt of the squeezed state in phase space is easy to achieve by letting it evolve, 
after its formation at tf, in the fixed, non-rotating trap. 

As a simple example let us consider the generation of squeezed vacuum states starting from the ground state of 
the initial trap, so that a = 0. To design the squeezing process we may follow a similar procedure as in the previous 
section, but minimizing the cost function 

F = 
b = 

which is minimal for h{tf) = 1 and h{tf) = 0 , so that b{tf) = 7 and b{tf) = 0 . 

Since, due to the centrifugal force during the rotation, the wave packet tends to spread first, the squeezed states 
with 7 > 1 may be achieved in shorter times than the ones needed without squeezing in the previous section. Figure 
[S] depicts the free parameters that optimize a rotation with a final squeezed state for the same parameters in the 
previous subsection, but 7 ^ = 3, and Fig. |4] the excess energy with respect to the target state. The excitation in 
a process with a final moderate squeezing is smaller than for the simple rotation without squeezing. Fig. [S] depicts 
the difference between the target value of the function b (proportional to the width of the wavepacket) and its actual 
value at final time for rotations without and with squeezing. Again, the minimizations change suddenly to a different 
solution that cannot satisfy the conditions at a critical time, see also Figs. [3]and|6l 


b{tfY 

6 - 7 -I-I, (40) 


VII. DISCUSSION 

We have worked out different schemes to perform fast rotations of a one-dimensional trap without any final excitation 
of the confined particle, which we have considered to be an ion throughout but could be a neutral particle as well by 
setting the proper trapping interaction. Apart from excitation-free rotations it is also possible to generate squeezed 
states in a controllable way. For an arbitrary trap, the fast processes could in principle be performed in an arbitrarily 
short time if an auxiliary harmonic potential with time dependent frequency could be implemented. In a simpler 
setting, where only the rotation speed may be controlled, the rotation time cannot be arbitrarily short, as demonstrated 
by inverse engineering or bang-bang approaches, and confirmed by optimal-control theory. Bang-bang and optimal 
control protocols provide useful information and time bounds but are difficult to implement experimentally due to the 
sudden kicks in the angular velocity of the trap. Smooth protocols designed by invariant-based inverse engineering 
have also been worked out. They achieve negligible excitations for times close to the minimum times given by optimal 
control theory. 

The analysis may be generalized for a two-dimensional trap but it becomes considerably more involved (2^ and will 
be considered separately. The ID approximation used here will be valid for total energies well below the transversal 
confinement energy A_l = fiw±. For the shortest final times considered in our simulations, excitation energies are 
never larger than 2huJo so that uj± ^ ujq would be enough for their validity. 

Rotations are elementary manipulations which together with transport, splitting, and expansions, may help to 
build a scalable quantum information architecture. In particular, they provide a mechanism for connecting sites by 
changing transport directions in 2D networks. Rotations have been demonstrated experimentally for trapped ions 
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[13 and improving the capability to control the parameters involved is feasible with state-of-the-art trapped-ion 
technology. To extend the present analysis to ion chains an approach similar to that in iilllll could be 
applied, working out the dynamical modes of the system and taking into account the dipole-dipole interaction due to 
the rotation of the charged particles. The present results set a first step towards accurately controlling rotating ion 
chains which would allow for fast reordering. 
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Appendix A: Wave functions 


The time-dependent wave functions evolving with the Hamiltonian ([8]) take the form 

(s|V;(t)) = ^c„e*“"«(s|<^4t)) (Al) 

n 

where the c„ are constant, 

{s\4'n{t)) = (A3) 

and $„(x) is the Hermite polynomial solution of the harmonic oscillator with angular frequency wq and energy 
eigenvalue (n -|- l/2)fiwo, ^nix) = 2 >r Note that p^4>„(s/6) is just a scaled state 

which corresponds to the n-th eigenstate of a trap with angular frequency uJojb'^. 
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